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ABSTRACT 

In this paper we present a temporal and spectral analysis of X-ray data from the XMM-Newton and Chandra 
observations of the ultrasoft and variable Seyfert galaxy RX J1301. 9+2747. In both observations the source clearly 
displays two distinct states in the X-ray band, a long quiescent state and a short flare (or eruptive) state which differs 
in count rates by a factor of 5-7. The transition from quiescent to flare state occurs in 1-2 ks. We have observed 
that the quiescent state spectrum is unpreccdcntcdly steep with a photon index F" ~ 7.1, and the spectrum of the 
flare state is flatter with V ~ 4.4. X-rays above 2 keV were not significantly detected in either state. In the quiescent 
state, the spectrum appears to be dominated by a black body component of temperature about ~30-40 eV, which 
is comparable to the expected maximum effective temperature from the inner accretion disk. The quiescent state 
however, requires an additional steep power-law, presumably arising from the Comptonization by transient heated 
electrons. Optical spectrum from the Sloan Digital Sky Survey shows Seyfert-like narrow lines for RX J1301. 9+2747, 
while the HST imaging reveals a central point source for the object. In order to precisely determine the hard X-ray 
component, future longer X-ray observations are required. This will help constrain the accretion disk model for RX 
J1301. 9+2747, and shed new light into the characteristics of the corona and accretion flows around black holes. 
Subject headings: accretion, accretion disks - - galaxies: active - - galaxies: individual (RX 
J1301.9+2747) — X-rays: galaxies 



1. INTRODUCTION 

Active galactic nuclei (AGNs) are thought to be pow- 
ered by supermassive black holes of Msh ~ 10 6 — 10 9 M@ 
accreting the surrounding gas (see [Rces 1984 f° r a re- 
view). They are also considered to be scaled- up versions 
of Galactic blac k hole binaries (BHBs, Mbh ~ 10 M@, 
iMcHardv et al.1 (J2006I ) and references therein). The 
rapid X-ray variability is one example of the similari- 
ties b etween these two types of systems fGierlinsk i et al.l 
120081 ). In Seyfert galaxies, variations of the X-ray con- 
tinuum emission over a t imescale from minutes to hours 
have been reported (e.g . lUlrich et alj|1997t IBoller et al.l 
Il997t [Ponti et al.l 120121 ). however, persistent giant and 
rapid variability appears to be fairly rare and its origin 
is still poorly understood. 

Soft X-ray excesses above an extrapolation of the un- 
derlying hard X-ray power-law is comm only observed in 
Type 1 AGNs and radio quiet quasars (jPiconcelli et al.l 
120051: IBianchiet aT1l2009D . The origin of this additional 
component is not clear, and may be the hig h-energy tail 
of th e AGN accretion disk emission (e.g. iGrupe et al.l 
119951 ). The problem using this explanation is that the 
temperatures of the soft X-ray excesses appear to fall 
within a narrow range (kT ~ 0.1-0.2 keV) from a sam- 
ple of AGNs containing a large range of black hole 
(BH) mass, which is difficult to ex plain using the stan- 
dard accretion disk m odels (e.g. lGicrlinski fc D one 2004; 
iCrummv et~aTl 120061) . 

So far there has been no convincing evidence for the 
presence of direct ac cretion disk emissi on seen in the X- 
ray spectra of AGNs. lYuan et all (J2010I ) reported a lumi- 
nous ultra-soft excess in the narrow line Seyfert 1 galaxy 
(NLS1) J1633+4718 from archival ROSAT spectra, and 
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found a lowest soft excess temperature of 32 eV among 
AGNs. This characteristic of the soft excess is likely an 
observational signature for the accretion disk emission. 
However, the blackbody nature of this emission needs to 
be tested further, utilizing higher quality X-ray data, as 
the ROSAT spectra (~ 0.1-2.4 keV) are less sensitive to 
constrain the harder powe r -law e mission above ^2 keV. 
Recently, iTerashima et al.l (|2012f) reported the discovery 
of a candidate 'ultrasoft' AGN, whose X-ray spectrum 
can be represented purely by a soft thermal component 
with a blackbody temperature of kT ^0.13-0.15 keV, 
by analog with the accretion disk dominated spectrum 
typically seen in the high/soft state of BHBs. Addi- 
tionally, the soft X-ray emission obtained shows spectral 
variability consistent with being caused by strong Comp- 
tonization. Interestingly, the obj ect was later o ptically 
confirmed to be a Type 2 AGN (|Ho et al.ll2012T) with a 
central BH mass as small as 1O 5 M . However, in their 
work they did not test in detail the possibility of the 
accretion disk emission as the origin for the soft excess. 
In this paper, we report results of new XMM- 
Newton and Chandra observations of RX 1301.9+2747 
at z=0.0237 (hereafter J1302), a highly variable 
and ultra-soft X-ray so urce in a post-starburst galaxy 
( Dcwaiig an et al.l [2000) . Our detailed analysis of the 
optical spectrum from the SDSS revealed that it is a 
Seyfert galaxy. The ultrasoft X-ray emission of J1302 
was confirmed in the new X-ray observations. In partic- 
ular, we found unusual giant flares in both XMM-Newton 
and Chandra light curves, accompanied by spectra hard- 
ening during the flare state. In the quiescent state, 
the spectrum appears to be dominated by a thermal 
blackbody component, whose temperature is compara- 
ble to the predicted maximum accretion disk tempera- 
ture. Throughout this paper, we assume a cosmology 
with H = 0.71, ftjw = 0.27, n A = 0.73. 

2. DATA ANALYSIS AND RESULT 



2.1. Optical Spectrum 

J1302 was spectroscopically observed by the SDSS in 
March 2007. Figure 1 shows the rest-frame spectrum for 
J 1302 (black line), which is dominated by the starlight 
of the host galaxy. To subtract the starlight and the 
nuclear c ontinuum, we follo wed the recipe described in 
detail in iDong et al.l ()2005t ). As seen in Figure 1, the 
galaxy starlight model (green line) gives a very good fit 
to the optical continuum (x 2 /dof = 3648/3208). After 
the subtraction of stellar absorption lines, we fitted the 
emission-line spectrum, represented by a blue line, by 
using Gaussians to derive the parameters of the emis- 
sion lines. [OIII] AA4959, 5007, [Nil] AA6548, 6583, Ha 
and [SII] AA6717, 6731 emission lines are clearly detected 
with S/N > 5, while H/3 line is only weakly distinguish- 
able with S/N ~ 1.4. The right panel in Fig. 1 displays 
the emission-line spectra, alongside the best fit Gaussian 
models. There is no apparent broad component of Ha 
line, and a narrow Gaussian with a line width of FWHM 
~ 240km s _1 can provide good fit. In order to verify the 
absence of the broad Ha line, we add an additional Gaus- 
sian to the narrow Ha, with the width fixed at 2000 km 
s _1 . The flux is allowed to vary in the fitting process. We 
found that the fitting was marginally improved by adding 
this component, and the S/N for the Ha broad compo- 
nent is only ~0.5. This suggests that the broad Ha line, 
if there is any, is extremely weak in J1302. The ratios of 
the narrow lines [OIII] A5007/H^ > 4.8 (using 3ct upper 
limit of H/3 line flux), and [Nil] A6583/Ha = 2.3, place 
J1302 into th e Sey fert regime on the BPT diagram of 
iKewlev et all (|2006f ). The flux ratios of [SII] /Ha and 
[OI] /Ha are 0.56 and 0.32, respectively, further strength- 
ening the Seyfert natu re of J1302 according to line ratio 
diagnostic diagrams of IKewlev et al.l (|2006l) . 

2.2. X-ray Observations 

J1302 was observed by XMM-Newton EPIC cameras 
in December 2000 with a total exposure time of 29 ks. 
It was detected ~7.3 arcmin away from the center of 
the field of view in the XMM-Newton imaging of the 
Coma cluster (ObsID 0124710801). The XMM-Newton 
data were reprocessed with the Science Analysis Software 
version 11.0.0, using the calibration files as of December 
2011. We used principally the PN data, which have much 
higher sensitivity, using the MOS data only to check for 
consistency. Spectra and light curves of source were 
extracted from a circular region with a radius of 40" 
centered at the source position for both PN and MOS 
cameras. Background spectra were made from source- 
free areas on the same chip using four circular regions 
identical to the source region. The epochs of high back- 
ground events were examined and excluded by using the 
light curves in the energy band above 12 keV. 

The Chandra pointing observation of J1302 was taken 
in June 2009 for about 5 ks. The data were processed 
with CIA0 (version 4.3) and CALDB (version 4.4.1), fol- 
lowing the standard criteria. Fig. 2 shows the Chandra 
X-ray contours of J1302, overlaid on the HST image in 
the B band. It is clear from the figure that the X-ray 
emission of J1302 is point-like. The center for the X-ray 
source is coincident with the optical nucleus of the galaxy 
(with a positional offset of ~0.1"). Given the subarcsec- 
ond spatial resolution of Chandra, we conclude that most 



of the X-ray emission, if not all, comes from the nuclear 
region of the galaxy, likely related to the AGN. 

2.3. X-ray light curve 

Light curves from the XMM-Newton and Chandra ob- 
servations are shown in Fig. 3. The source exhibits large- 
amplitude count rate variations in both observations. It 
can be seen from the PN light curve that there is a giant 
flare with count rates rising by a factor of 5 times the 
average value, having a duration of ~2 ks. The X-ray 
flux then declines to a relatively steady state. Such X- 
ray flare is confirmed by the MOS light curves, in which 
a possible decline of another flare is also recorded at 
the beginning of the MOS observation. The time inter- 
val of the two flares is about 17 ks. Interestingly, similar 
flare is seen in the Chandra light curve, with count rates 
increasing by a factor of 7 within ~1 ks. The similar 
amplitudes of flares in the XMM-Newton and Chandra 
observations which span ~ 9 years suggest that the flare 
behaviour in J1302 seems persistent on time scale of ~ 
decade. The spectral variability during the flare will be 
explored in detail in Section [ 



2.4. X-ray Spectra 

As both XMM-Newton and Chandra observations show 
peculiar flare behaviours, we attempt to quantify the 
spectral variability during flares by dividing the data into 
high and low flux intervals, using count rate thresholds 
of 0.35 counts s" 1 for XMM -PN and 0.08 counts s" 1 for 
Chandra, respectively. For simplicity, we classify the 
data above the count rate thresholds as belonging to the 
flare state, and that which falls below, to the quiescent 
state. The spectra data were grouped in the following 
manner: data from the XMM-Newton had at least 20 
counts per bin to ensure the x 2 statistics, and the Chan- 
dra data had at least 3 counts per bin and utilized the C- 
statistics which was adopted for minimization. Spectral 
fitting was performed using the XSPEC (Version 12.6) and 
limited to the 0.2-2 keV range for XMM-Newton, since 
the emission is background dominated above that energy 
range. The Chandra data was fitted in the energy range 
between 0.3 and 3 keV. Throughout the model fittings, 
the Galactic column densi ty was considered and fixed at 
jyGai = 0.75 x iq20 cm -2 (jKalberla et al.|[200l . 

Spectral variability is clearly present between the two 
states with the source being much harder in the flare 
state. To illustrate the differences in the spectral slope 
between the two states, we show the spectra together 
with an absorbed power-law model in Figure 4. While 
the fit is generally acceptable for the two states, as con- 
firmed by the reduced \ 2 value (see Table 1), the photon 
indices obtained from a power-law fit are, however, ex- 
tremely steep with T = 4.4±°j for N H = 4.3t 2 ;§ x 10 20 



cm" 2 , and T = 7.1Io:? for n h = 3.61^ x 10 2U cm" 2 , 
for the XMM flare and quiescent state, respectively. 
The photon indices belong to the steepest values ob- 
tained from AGN X-ray spectra. For comparison, the 
mean photon index in the 0.2-2.0 keV band is ~ 2.9 
for a sample of soft X-ray selected AGN s observed with 
the ROSAT RASS flGrupe et all l20"Toft . The result 
of this comparison indicates that the X-ray spectrum 
of J1302 is extremely soft compared to other AGNs. 
The absorption-corrected luminosity in the 0.5-2.0 keV 



range for this simple power- law model is 6.7 x 10 and 
2.8 x 10 40 erg s" 1 for the XMM flare and quiescent 
state, respectively. 

In order to further investigate the spectral variability in 
J1302, we then attempted to fit the spectra with a black- 
body (bbody in XSPEC) or Multiple Color Disk model 
(MCD, diskbb in XSPEC), an d a thermally Com ptonized 
disk model (compTT in XSPEC. iTitarchuk lfl994l) . both of 
these alternative models ha ve been used to fit the spec- 
tra of Galactic BHBs (e.g.. iDone et al.l 12007ft. and the 
soft X-ray excess emission in AGNs (e.g.. iPorquet et aTl 
120041: IPatricket aTll2012D . The diskbb model inte- 
grates over the surface of accretion disk to form a mul- 
ticolor blackbody spectrum, and compTT is an analytic 
model that self-consistently calculates the spectrum pro- 
duced by the Comptonization of soft seed photons in 
a hot corona above the accretion disk. The physical 
parameters of the compTT model are: the soft photon 
temperature (fcTo), the temperature of the Comptoniz- 
ing electrons (kT e ), the plasma scattering optical depth 
(r). For our fitting with the compTT model, a disk geom- 
etry was assumed for the comptonizing region, and the 
seed photons were assumed to follow Wien's law with 
a temperature of 22 eV (the expected disk temperature 
in section I3.2p . Because the temperature and optical 
depth of the Comptonizing plasma are strongly coupled 
(both are equally involved in shaping the spectrum) and 
thus cannot be constrained simultaneously, we fixed the 
plasma temperature at 20 keV and obtained constraints 
on the optical deptlfl. The single Comptonized model 
yields consistent fitting results with the previous simple 
power-law model for the spectrum at both states. The 
Compton optical depth is r = O.I6+0Q5 and r < 0.03 
for the flare and quiescent state, respectively. For the 
spectrum in the XMM-Newton flare state, a multicolor- 
disk blackbody gives equivalent fit, which is statistically 
better than the simple blackbody emission. 

The spectrum at the XMM-Newton quiescent state, 
however, shows an excess of emission at energies above 
~ 0.7 keV when fitted with a thermal model (bbody or 
diskbb). The addition of a power-law to the model im- 
proves the fit with very high statistical significance (x 2 
decreased by 13.9 for two extra parameters, at a 99.98% 
level according to _F-tcst). The power- law component 
contributes ^15% of the total luminosity in the 0.3-2 
keV band. In this case, we obtain an effective black- 
body temperature of /cTbb = 43lg eV, comparable to 
the seed photon temperature assumed in the compTT 
model. Although with large uncertainties, the addi- 
tional power-law component is relatively steep with pho- 
ton index T = 4.3^'g, and it is close to what is observed 
in the flare state. The spectral fitting results for the PN 
data, alongside the observed flux and intrinsic luminos- 
ity in the 0.5-2 keV range for each model, are shown in 
Table 1. Note that we used the same models to fit the 
MOS data and found that the results agree well with the 
PN data. 

The Chandra spectra were fitted with the same mod- 
els used in the XMM-Newton observation and the results 

2 Leaving the plasma temperature as a free parameter yields 
kT e = 21(< 27) keV and r < 1.3 for the XMM quiescent state 
spectrum, but both parameters cannot be constrained by the data 
during the flare. 



are listed in Table 1. During the first run we found that 
the photon indices derived from the simplest power-law 
model are systematically flatter than the values for the 
XMM-Newton data. The spectrum during the flare can 
be well fitted by a power-law model (C/dof=22.5/22). 
On the other hand, a power-law model is not sufficient 
to fit the data at the Chandra quiescent state. The ad- 
dition of a soft thermal component with respect to the 
power-law model improves the fit significantly (C value 
decreases by ^9 for two extra parameters, corresponding 
to a significance level of 98.7%). The resulting best-fit 
photon index is flatter, with T = 3.5^"'q. Similarly, we 
obtain an effective disk temperature, when fitted with a 
blackbody, of kT BB = 29±^ eV. This value is slightly 
lower than the blackbody temperature for the XMM- 
Newton quiescent state, but the parameter is loosely con- 
strained due to the poor statistics of the Chandra data. 
The power-law component in this case contributes ~ 40% 
of the total luminosity in the 0.3-2 keV, indicating a pos- 
sible change of the power-law or the blackbody emission 
between Chandra and XMM-Newton observations. 

Although a Comptonized model as opposed to a power- 
law model also provides a good fit for the Chandra flare 
state spectrum, it is not sufficient to fit the data at the 
quiescent state. The addition of an extra hard power-law 
component is needed at a significance level of 98.7% ac- 
cording to F-tcst. The best-fitted optical depth r for the 
compTT model is not significantly different from that ob- 
tained with the XMM-Newton data. The unabsorbed 
luminosity in the 0.5-2 keV band, based on the best-fit 
power-law model and the diskbb+power-law model for 
the Chandra flare and quiescent state, is 5.1 x 10 41 and 
4.4 x 10 40 erg s _1 , respectively. 

3. DISCUSSION 

3.1. AGN characteristics in RX 31301. 9+2747 

X-ray observation with Chandra , which has superb 
spatial resolution of ~ 0.5", revealed the presence of an 
AGN in J1302: the center of the bright unresolved X-ray 
emission coinciding with the optical point-like nucleus 
of the galaxy (with a position offset of ~ 0.1"). Both 
the Chandra and XMM-Newton observations show that 
it has Seyfert-like X-ray luminosity of ~ 10 41 erg s _1 in 



the energy range of 0.5-2 keV ( ~ 10 erg s 



the 



0.2-2 keV band ), and rapid X-ray variability down to 
a time-scale ~ 1 ks. Additionally, the optical spectrum 
of J1302 displays Seyfert-like narrow emission line ratios. 
All these observational facts point to the presence of an 
AGN (Seyfert nucleus) in this galaxy. This can be fur- 
ther supported by the point-like appearance of a nuclear 
source fro m the HST imaging o bservations with ~ 0.1" 
resolution (|Caldwell et all 119991 ) . 

Bas ed on ROSAT X-ray observations. IDewangan et al.l 
(|2000l) argued for the presence of an AGN in J1302, a 
view that is consistent with ours on the basis of the new 
observations. However, based on their observed optical 
spectrum, Dewangan et al. conclude that the galaxy nu- 
cleus is more like a LINER. In this paper, we have care- 
fully modeled the host galaxy's starlight, especially the 
stellar absorption features, and subtracted them from the 
new SDSS spectrum, enabling us to accurat ely measure 
the w eak AGN emission lines in J1302 (e.g.. lDong et al.l 
120051) . thus more resolutely confirming the Seyfert na- 



ture based on the line ratio diagnostics. 

The lack of detectable broad permitted lines prevents 
us from estimating the central BH mass of J1302 using 
conventional lincwidth-luminosity-mass scaling relation. 
In the SDSS spectral fitting, we found the strongest nar- 
row line, [OIIIJA5007, is marginally resolved with Gaus- 
sian a = 58 ± 9 km s _1 (after correcting for the instru- 
ment resolution). Using the width of the [OIIIJA5007 line 
as a proxy for the stellar velocity dispersion of the host 
galaxy, we obtained a BH mass of M bh = 8 x 10 5 M^ w ith 
an intrinsic scatter of 0.5 dex (e.g. JXiao et al.ll2011[) . 

The bolometric luminosity for J1302 can be estimated 
from the optical continuum luminosity. We retrieved 
the high resolution HST/WFPC2 images of J1302 in 
the B (F450W filter) and I (F814W filter) passbands 
from the HST archive. The HST observations (datasct 
U39D0301M-U39D0304M) were made in July 1997, with 
two 600 s exposures in the B band and two 400 s ex- 
posures in the I band, respectively. The images were 
processed usingthc standard HST pipeline routines in 
IRAF/STSDAS_j. We then performed two-dimensional 
profile decompositions of this galaxy with the code 
GALFIT (version 3.0, Peng et al. 2010). Our model con- 
sists of an exponential disk component, a bulge compo- 
nent, and an unresolved central point source for the nu- 
clear AGN emission. In our GALFIT modeling, the point- 
spread function w as generated by the Tiny Tim software 
(|Krist et al.lll995l) . Note that the bulge for this galaxy 
displays a box/peanut shape, which is commonly seen 
in edge-on barred disk galaxy, consequently we added a 
boxiness parameter to the bulge profile in GALFIT. The 
model generally matches the data well (xl ~ 1.02). The 
inferred flux for the central point source in the B-heaid 
is Mb = -15.8, corresponding to a nuclear luminosity of 
vL v b = 8 x 10 41 erg s _1 . For comparison, the B-band 
bulge luminosity of this galaxy derived from the GALFIT 
decomposition is Lb, bulge ~6x 10 42 erg s _1 . If indeed 
the nuclear emission comes from AGING, we estimate the 
bolometric luminosity to be 1 x 10 43 erg s _1 using the 
i?-band luminosity for the central point source by adopt- 
ing a bolometric correction of 13 (|Marconi et al.ll2004D . 
For a BH mass of 8 x 10 5 M Q , the accretion rate in Ed- 
dington unit is L/Ledd ~ 0.1, which suggests J1302 is 
accreting at high Eddington ratio. 

3.2. Extremely Soft X-ray Emission 

One of the remarkable features of J1302 is the extreme 
softness of the X-ray spectra. The best-fitted power- 
law index for the spectrum in the XMM-Newton qui- 
escent state (r ~ 7) is one of the steepest soft X-ray 
photon indices am ong AGNs (e.g., iGrupe et al.l 119951 : 
iBoller et al.ll201ll ). Understanding the origin of the ul- 
trasoft X-ray emission will help to pin down the nature of 
the s ource. Some AGN s such as NLSls can be very soft 
(e.g.. IBoller et al.lll996l:lMiddleton et alj|2007l) . showing 
a strong soft X-ray excess over an underlying power-law 
component. The strength of soft excess can be quanti- 
tatively described as the ratio of flux at 0.5 keV to the 
power law ext rapolation of the fittin g to the spectrum 
above 2 keV (jMiddleton et al.ll2007l) . Since no signifi- 

3 http://www.stsci.odu/institutc/softwaro_hardwarc/stsdas 

4 Note that we have not accounted for the dust extinction in 
estimating the luminosity. 



cant hard X-ray emission above ~2 keV was detected 
in the XMM observation of J1302, we estimated 90% 
confidence upper limit on the count rates in 2-7 keV, 
~ 6.7 x 10~ 4 counts s _1 , and converted it to an upper 
limit on the extrapolated flux at 0.5 keV. XSPEC simula- 
tions of models with power-law T = 2, 2.5 and 3 show 
that the lower limits on the ratios defined above are 37, 
15 and 6.1 , respectively. Note tha t the ratios for a sample 
of NLSls (JMiddleton et al.ll2007l) are usually found to be 
less than 10. Thus, the soft emission relative to that at 
the hard X-rays in J 1302 is extremely strong compared 
to other AGNs. 

The origin of soft excess is still unclear. The spec- 
trum at the XMM quiescent state can be fitted well 
by a Comptonizcd model and a blackbody plus power- 
law emission equally. Interestingly, the fitted black- 
body temperature (fcTeB = 43^ eV) is much lower 
than t he canonical values o f ~ 0.1-0.2 keV found for 
AGNs flCrummv et al.l l2006f) . Sta ndard accretion disk 
models (jShakura fc Sunvaevl H973I ) give a maximum ef- 
fective temperature of the accreting material hT^^x ~ 
ll^rh/Mg) 1 / 4 eV, where to is mass accretion rate in Ed- 
dington unit and M 8 = M B h/10 8 M q . Using Mbh = 8 x 
1O 5 M and to ~ 0.1 for J1302, we obtained fcT max ~22 
eV, which is comparable to the fitted blackbody tem- 
perature. Note that though with larger uncertainties, 
the lower fitted temperature for the Chandra quiescent 
state data (&Tbb = 29l^g eV) is more compatible with 
the predicted maximum disk temperature. Therefore, 
the ultrasoft X-ray emission in J1302 may be connected 
with the direct thermal emission from the accretion disk. 
Another constraint on the X-ray spectra due to thermal 
disk emission can come from the optical/UV data. The 
optical B-band flux of nuclear point source from the HST 
observation is, however, much higher than the extrapo- 
lated MCD flux in the optical (Mb ~ —12.7), about one 
order of magnitude difference. The difference cannot be 
explained by the contamination from nuclear star clus- 
ters, as they have ty pical absolute J-ba nd magnitudes be- 
tween — 14 and —10 (|B5ker e t al. 2002]), much lower than 
the observed J-band flux of J1302 nucleus (Mi ~ —16.8, 
obtained with the same GALFIT decomposition of the 
HST/WFPC2 image as detailed in Section EH). To fur- 
ther investigate this difference we need a more realistic 
disk spectral modeling, and fit to a broader band data, 
which is beyond the scope of this paper. 

3.3. Unusual X-ray variability 

Another unusual feature of J 1302 is that it clearly 
shows two distinct states in the X-ray, a flare (or erup- 
tive) state and a quiescent state. The amplitudes of the 
flare in the Chandra and XMM-Newton light curves look 
very similar, with the count rates increased by a factor of 
5-7 within ~l-2 ks. Both the XMM-Newton observation 
in 2000 and the Chandra observation in 2009 detected 
rapid flares, suggesting that the flare itself appears repet- 
itive and occurs very frequently in the object. In fact, 
ROSAT PSPC observations also found that the object is 
highly variable and demon strates a rapid flare ey ent in 
light curve lasting ~1.3 ks (|Dewangan et al. 2000). The 
rapid energetic flare in J 1302 is fairly rare among Seyfcrt 
galaxies and quasars, although some extreme variations 
have also been found in objects such as IRAS 13224- 



3809 ( [Boiler et alJH997h and PKS 05 58-504 dWang et al.1 
120011) . As discussed by IWang et~all (|2001[ ). a magnetic 
coronal model in which electrons in the corona are con- 
tinuously heated by magnetic reconnection, can produce 
rapid energetic X-ray flare. A realistic physical model to 
explain a rapid energetic flare in AGNs, and the compar- 
ison with Galactic BHBs are beyond the scope of this 
work and will be presented elsewhere. 

The spectral variability is clearly present between the 
two flux states with the source being much harder in 
the flare state. This behavior is markedly contrast 
to what is commonly seen in AGNs and BH Bs (e.g., 
iMarkowitz fc Edelson I [20041; iDone et all 120071) . Clearly 
the spectral variability cannot be explained by any simple 
spectral model. One possibility is that the spectrum in 
the quiescent state is represented by a relatively stable 
thermal emission from accretion disk, and the spectral 
variability is caused by the flux variations of a second ad- 
ditional component such as strongly Comptonized power- 
law emission. This perhaps explains why the photon in- 
dex for the power-law emission in the XMM-Newton flare 
state (which dominates the spectrum) is consistent with 
the additional power-law component in the quiescent 



state, whose strength is relatively weak. In this case, 
the underlying power-law emission is still steep (r ^4) 
compared to other AGNs, the nature of which remains 
understood. Note that similar spectral change has also 
been seen in the Chandra observation of J1302, though 
the spectra statistics for the data is low. Longer X- 
ray observations are required to examine the presence of 
persistent flares and to investigate the origin of spectral 
variability during flare periods. This can in turn shed 
new lights on the characteristics of corona and accretion 
flows around BHs. 
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TABLE 1 
Spectral fitting results for the XMM-Newton and Chandra observation at different states 
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a Spectral model (as given in XSPEC) multiplied by a neutral absorption (wabs) with column density Nh- 

photo-electric absorption model using Wisconsin-Morrison & McCammon (1983) cross-sections. 
b Plasma optical depth in the CompTT model. 



wabs is the 



energy range of 0.5—2 keV, in units of 10 4t) erg s _1 . 
d The column density was fixed to Galactic value N^| al if the fitting yields a Nh <N^ al . 
e The spectrum of the seed photons is assumed to be Wien law with a temperature of 22 eV (see Section 13.21 for details). We 

fixed the plasma temperature at 20 keV and obtained constraints on the optical depth. 
f Pegged at the minimum value allowed in XSPEC. 
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Fig. 1. — Illustration of the continuum and emission-line fittings of the SDSS spectrum. Left panel: observed spectrum (black), stellar 
continuum model (green) and residual (blue) which is used to fit the emission lines. Right panel: a zoomed-in view of the emission-line 
profile fitting for the Hq+[NII] and the [SII] doublet lines. The inset in the left panel shows a zoomed-in view of the H/3 and [OIII] region. 
Gaussian line models are plotted in cyan line, and the residual is shown in the lower panel. In the emission-line fits, the line ratios of 
[OIII] and [Nil] doublet lines are fixed to the theoretical values. 
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Fig. 2. — Contours of the Chandra image (green) overlaid with the HST S-band image. The direction of north is up and east is left. 
White plus marks the posit ion of the central p oint source from the HST imaging. For comparison, the center of the radio emission obtained 
with the Very Large Array (Miller ct al. 2009J) is shown in cyan cross. 
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Fig. 3. — XMM-Newton and Chandra light curves for J1302. Top panel: XMM-PN and the summed MOS1+MOS2 light curves, with 
a time bin size of 100 s. The PN and MOS background light curves are also shown for comparison (lower panel). Bottom panel: Chandra 
light curve with the same bin size as XMM-Newton. The dotted lines represent the count rate thresholds dividing the data into flare and 
quiescent state, which are 0.35 counts s — 1 for XMM-PN, and 0.08 counts s — 1 for Chandra , respectively. 
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Fig. 4. — Top panel: Spectra of the XMM-Newton flare state (black) and the quiescent state (red) for J1302. Only the PN data are 
shown for clarity. The solid lines are the simple power-law model fits for both states, and the corresponding residuals are shown in the 
lower panels. Bottom panel: as top panel, but showing for the Chandra data. 



